ADA020471 


TEXAS  A&M  UNIVERSITY 


COLLEGE  STATION,  TEXAS  77843 


m MM 


TEXAS  A&M  UNIVERSITY 
PROJECT  THEMIS 


Technical  Report  No.  50 

STATISTICAL  PERT:  DECOMPOSING 

A PROJECT  NETWORK 


by 


R.  L.  Sielken  Jr.  and  N.  E.  Fisher 


ATTACHMENT  I 


STATISTICAL  PERT:  DECOMPOSING  A PROJECT  NETWORK 


by 


R.  L.  Sielken  Jr.  and  N.  E.  Fisher 


THEMIS  OPTIMIZATION  RESEARCH  PROGRAM 
Technical  Report  No.  50 


INSTITUTE  OF  STATISTICS 
.Texas  ASM  University 
November  1975 


Research  conducted  through  the 
Texas  ASM  Research  Foundation 
and  sponsored  by  the 
Office  of  Naval  Research 
Contract  N00014-68-A-0140 
Project  NR047-700 


Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of 
the  United  States  Government. 


This  document  has  been 
approved  for  public  release 
and  sale;  its  distribution 
is  unlimited. 


ATTACHMENT  II 


TABLE  OF  CONTENTS 


Program 

BREAKUP 

General  Description 
Specific  Input  Instructions 
Dimension  Restrictions 
Sample  Network 
Sample  Input 
Sample  Output 

Sample  Problem:  Schematic  Representation 

Program  Listing 
LOOP 

Specific  Input  Instructions 
Dimension  Restrictions 
Sample  Problem 
Program  Listing 


L_ 


ABSTRACT 


Statistical  PERT  is  a new  procedure  for  obtaining  information 
about  the  distribution  of  a project's  completion  time  when  the 
project  is  comprised  of  a large  number  of  activities  and  the  time 
required  to  complete  an  individual  activity  once  it  can  be  begun 
is  a random  variable.  The  project  is  represented  as  an  acyclic 
network  whose  arcs  correspond  to  the  project  activities.  This 
network  is  simplified  by  replacing  various  activity  configurations 
by  single  equivalent  activities  and  then  decomposed  into  several 
subnetworks.  — <Hie  distribution  and  momenta  of  each  subnetwork's 
completion  time  ane  bounded  and  approximated  on  the  basis  of  two 
points  from  each  .activity's  completion  time  distribution  by  using 
some  mathematical  programming  techniques  and  a new  result  in  the 
theory  or networks.  The  project's  completion  time  distribution  is 
then  approximated  by  combining  the  approximate  subnetwork 
distributions. 

^ This  report  documents  two  computer  programs.  The  first  program 
BREAKUP  decomposes  a project  network  into  several  subnetworks  which 
are  connected  in  either  series  or  parallel  in  the  project  network. 

The  second  program  LOOP  checks  a given  project  network  for  loops 
(cycles)  since  any  loop  would  contradict  the  assumed  acyclic  structure 


of  the  project  network. 
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BREAKUP 


General  Description; 

This  program  "breaks  up"  a network  into  a set  of  subnetworks 
which  can  be  linked  together  either  in  series  or  parallel  to  yield  the 
given  network.  This  breakup  is  complete  in  the  sense  that  none  of  the 
subnetworks  in  the  set  can  be  further  broken  up. 

The  basic  breakup  procedure  involves  two  main  subroutines,  BUNDLE 
and  CUT.  BUNDLE  partitions  the  activities  in  a given  network  or 
subnetwork  into  parallel  subnetworks  connecting  the  network's  source 
and  sink.  CUT  identifies  the  cut  nodes  in  a given  network  or 
subnetwork  and  then  identifies  the  sets  of  activities  between  each  of 
the  consecutive  cut  nodes.  This  series  of  activity  sets  represents 
a breakup  of  the  given  network  or  subnetwork  into  subnetworks  in 
series.  The  complete  breakup  of  the  given  network  is  the  following 
sequential  procedure: 

(1)  Use  BUNDLE  to  identify  the  parallel  subnetworks  connecting 
the  source  and  sink  of  the  network.. 

(2)  Use  CUT  separately  on  each  parallel  subnetwork  identified 
in  the  previous  step  - (1)  or  (3)  and  breakup  the 
parallel  subnetwork  into  subnetworks  in  series.  If  no 
such  breakup  is  possible  for  a parallel  subnetwork,  that 
subnetwork  is  not  considered  again.  If  no  new  series 
subnetwork  is  identified  in  this  entire  step,  stop. 

Otherwise  go  to  (3). 

(3)  Use  BUNDLE  separately  on  each  series  subnetwork  identified 
in  (2).  If  BUNDLE  cannot  breakup  a series  subnetwork, 
that  subnetwork  is  not  considered  again.  If  no  new 
parallel  subnetwork  is  identified  in  this  entire  step. 


A schematic  example  of  this  procedure  accompanies  the  sample  problem. 


/ 
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Specific  Input  Instructions: 


Card  1.  Col.  1-3: 
Col.  4-6: 


Col.  7-9: 


Col.  10-12: 


The  number  of  arcs  in  the  network,  Format  (13). 

The  number  of  the  node  which  is  the  source  node. 
Format  (13) . 

The  number  of  the  node  which  is  the  sink  node. 
Format  (13). 

The  largest  node  number  in  the  network.  Format  (13) 


For  each  activity  one  card  with: 

Col.  1-3:  The  activity's  number.  Format  (13). 

Col.  4-6:  The  activity's  origin  node  number,  Format  (13). 

Col.  7-9:  The  activity's  terminal  node  number,  Format  (13). 

The  activities  and  nodes  may  be  numbered  in  any  way  and  may  be  read  in 
in  any  order. 

Dimension  Restrictions: 

Currently  this  program  will  accomodate  a network  that  has  a 
maximum  of  300  nodes  and  a maximum  of  300  arcs.  It  can  store  a 
maximum  of  100  subnetworks  with  each  subnetwork  having  a maximum 
of  300  arcs. 

This  program  is  written  in  FORTRAN  G. 
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| STAGE  1 BREAKUP 

SUBNETWORK  2 IS  COMPOSED  OF  SUBNETWORKS: 

7,  8,  9,  10, 

IN  SERIES 

f 

SUBNETWORK  3 IS  COMPOSED  OF  SUBNETWORKS: 

11,  12,  13, 

IN  SERIES 

SUBNETWORK  4 IS  A MINIMUM  NETWORK 

I IT  IS  COMPOSED  OF: 

SOURCE  NODE  * 1 

SINK  NOOE  * 18 

ARC  SI  ARC  I TIARC) 

29  1 18 

SUBNETWORK  5 IS  COMPOSED  OF  SUBNETWORKS: 

14,  15,  16, 

| 

IN  SERIES 

I 

I 

SUBNETWORK  6 IS  A MINIMUM  NETWORK 
IT  IS  COMPOSEO  OF: 

(SOURCE  NOOE  ■ 1 

SINK  NODE  * 18 

ARC  SI  ARC)  TIARC) 

1 1 18 

1 


STAGE 


2 BREAKUP 


7 


SUBNETWORK  7 IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  » 1 

SINK  NODE  * 7 

ARC  S(ARC)  T( ARC ) 

2 1 7 


i 


SUBNETWORK  8 IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  « 7 

SINK  NOOE  * II 

ARC  S( ARC)  T ( ARC ) 


3 7 8 

A 7 9 

6 8 10 

8 8 11 

5 9 8 

7 9 10 

9 10  11 


| 

T 

fe 

SUBNETWORK  9 IS  COMPOSEO  OF  SUBNETWORKS: 

X 

I 

17,  18, 

IN  PARALLEL 

■i  I 


SUBNETWORK  10  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  > 17 

SINK  NOOE  * 18 

ARC  SI  ARC)  T( ARC ) 

19  17  18 


SUBNETWORK  11  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSEO  OF: 

SOURCE  NODE  > 1 


SINK  NOOE 


A 


8 


S ( ARC)  T( ARC ) 
l 2 

1 3 

2 4 

3 2 

3 4 

I 

SUBNETWORK  12  IS  COMPOSED  OF  SUBNETWORKS: 

19.  20. 

f 

IN  PARALLEL 

i 

i 

i 

SUBNETWORK  13  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSEO  OF: 

SOURCE  NODE  « 6 

SINK  NODE  * 18 

ARC  SI  ARC)  TIARC) 

28  6 18 


SUBNETWORK  14  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  ■ 1 

SINK  NODE  « 19 

ARC  SI  ARC)  TIARC) 

30  l 19 


SUBNETWORK  IS  IS  COMPOSED  OF  SUBNETWORKS: 


j 21.  22.  23. 


IN  PARALLEL 


SUBNETWORK  16  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSEO  OF: 

SOURCE  NODE  > 20 

SINK  NODE  * 18 

ARC  SI  ARC)  TIARC) 

34  20  18 
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STAGE  3 BREAKUP 

I 

SUBNETWURK  17  IS  CCMPOSED  OF  SUBNETWORKS: 
24,  2b» 

IN  SERIES 

SUBNETWORK  18  IS  COMPOSED  OF  SUBNETWORKS: 
2t>»  27,  28, 

IN  SERIES 


9 


I 


SUBNETWORK  19  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  * 4 

SINK  NODE  * 6 

ARC  SIARC)  T( ARC  I 
27  4 6 

SUBNETWORK  20  IS  COMPOSED  OF  SUBNETWORKS: 

29,  30, 

IN  SERIES 

SUBNETWORK  21  ISA  MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  « 19 

SINK  NODE  * 20 

| 

ARC  SIARC)  TIARC) 

31  19  20 

SUBNETWORK  22  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  > 19 

SINK  NODE  * 20 

ARC  SIARC)  TIARC) 

32  19  20 


I 


SUBNETWORK  23  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF S 


10 


SOURCE  NODE 

« 19 

SINK 

NODE 

* 20 

ARC 

S( ARC! 

TIARC) 

33 

19 

20 

4 BREAKUP 


SUBNETWORK  24  I S A MINIMUM  NETWORK 
IT  IS  COMPOSED  OFS 

SOURCE  NODE  • 11 

SINK  NODE  * 12 

ARC  SI  ARC ) TIARC) 

10  11  12 

SUBNETWORK  25  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  » 12 

SINK  NODE  * 17 

ARC  SI  ARC  I T I ARC  ) 

18  12  17 


SUBNETWORK  26  I S A MINIMUM  NETWORK 
IT  IS  COMPOSED  OF: 

SOURCE  NODE  « 11 

SINK  NODE  * 13 

ARC  SI  ARC)  TIARC) 

11  11  13 


SUBNETWORK  27  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSEO  OF: 


SOURCE  NOOE 

» 13 

SINK 

NODE 

» 16 

ARC 

SI  ARC ) 

TIARC) 

12 

13 

14 

13 

13 

15 

16 

14 

16 

15 

15 

16 

14 

15 

14 

SUBNETWORK  28  IS  A MINIMUM  NETWORK 
IT  IS  COMPOSEO  OF: 


SOURCE  NOOE  * 16 


1 


SINK  NODE  = 17 

ARC  S(ARC)  TIARC) 

17  16  17 

SUBNETWORK  29  IS  A MINIMUM  NETWORK 

IT  IS  COMPOSED  OF: 

SOURCE  NODE  = 4 

SINK  NODE  * 5 

ARC  S( ARC ) T ( ARC ) 

25  4 5 

SUBNETWORK  30  IS  A MINIMUM  NETWORK 

IT  IS  COMPOSED  OF: 

SOURCE  NODE  * 5 

SINK  NODE  * 6 

ARC  SI  ARC!  TIARCI 

26  5 6 
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SAMPLE  PROBLEM:  SCHEMATIC  REPRESENTATION 

BREAKUP 

Arcs  are  designated  as  follows: 

© KD 

The  arc  number  is  written  above  or  to  the  left  of  the  arc  activity  line. 
The  arrowhead  indicates  the  direction  of  the  activity.  The  circled 
numbers  are  node  numbers.  In  this  case  the  arc  depicted  is  arc 
number  10,  flowing  from  left  to  rJght  with  beginning  node  5 and 
terminating  node  7. 

Subnetworks  are  designated  as  follows: 


The  subnetwork  number  is  written  above  or  to  the  left  of  the  network 
activity  direction  line.  The  arrowhead  Indicates  the  direction  of  the 
subnetwork  activity.  The  numbers  enclosed  in  ovals  are  source  and 
sink  nodes.  In  this  case  the  subnetwork  depicted  is  subnetwork  11, 
flowing  from  left  to  right  with  source  node  12  and  sink  node  15.  A 
solid  activity  direction  line  indicates  the  subnetwork  is  a minimum 
network.  A broken  line  indicates  the  subnetwork  must  be  considered 
by  at  least  one  more  subroutine  of  the  BREAKUP  program. 

CD ><D 

The  above  figure  indicates  subnetwork  8,  flowing  from  left  to  right 
with  source  node  6 and  sink  node  17,  is  not  yet  a minimum  network. 

In  the  BREAKUP  diagrams,  the  decomposition  of  the  initial  network 
is  finished  when  all  subnetworks  are  minimum  networks,  e.g.  all  lines 


are  solid. 
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PROGRAM  BREAKUP 


BREAKUP  IDENTIFIES  PARALLEL  SUBNETWORKS  AND  SERIES  SUBNETWORKS 

MORE  THAN  ONE  NETWORK  MAY  BE  DECOMPOSED  DURING  A RUN. 

SEPARATE  THE  DATA  FOR  EACH  SUBNETWORK  BY  A BLANK  CARD 


THE  FOLLOWING  IS  AN  ALPHABETICAL  LISTING  UF  THE  VARIABLES  AND 
ARRAYS  THAT  ARE  USED  IN  THIS  MAIN  PROGRAM  AND  ITS  SUBROUTINES 


ARCS  * THE  NUMBER  OF  ARCS  IN  THE  SUBNETWORK 
BDNUM(I)  * THE  BUNDLE  NUMBER  TO  WHICH  NODE  I IS  ASSIGNED 
CHECK  * ARRAY  USED  TO  STORE  ARCS  HAVING  THE  SOURCE  AND  THE 
SINK  AS  THEIR  ONLY  NODES 

CTNSUB  « THE  TEMPORARY  NUMBER  OF  SUBNETWORKS  FOUND  IN  THE 
PREVIOUS  STFP 

LNODEN  = THE  LARGEST  NODE  NUMBFR  BEING  READ  IN 
MAXND  * THE  LARGEST  NODE  NUMBER  THAT  HAS  ALREADY  BEEN  ASSIGNED 
AT  LEAST  TEMPORARILY  TO  A BUNDLE 
N ARCS S(  I ) « THE  NUMBER  OF  ARCS  IN  SUBNETWORK  I 
NSUB  « THE  TOTAL  OF  SUBNETWORKS  THUS  FAR 
NTARC  = THE  NUMBER  CF  ARCS  IN  THIS  SUBNETWORK 
NUMBD  * THE  NUMBER  OF  BUNDLES  CREATED 
SID  * THE  STARTING  NGOE  FOR  ARC  I 
SINK  » THE  NOOE  NUMBER  CORRESPONDING  TO  THE  SINK 
SINKS! I ) - THE  SINK  IN  SUBNETWORK  I 
SnURC(I)  * SOURCE  NODE  IN  SUBNETWORK  I 
SOURCE  * THE  NOOE  NUMBER  CORRESPONDI MG  TO  THE  SOURCE 
STEP  * STAGE  NUMBER 

SUBNET! I» J)  « THE  ITH  ARC  IN  THE  JTH  SUBNETWORK 
SUMARC  * THE  CURRENT  NUMBER  OF  ARCS  IN  SUBNETWORK  NSUB 
T(I)  « THE  TERMINATING  NODE  FOR  ARC  I 

TARC  = THE  SUBNETWORK  WITHOUT  THE  ARCS  INVOLVING  NODE  K 
TLNDEN  « TEMPORARY  LARGEST  NODE  NUMBER 

TNSUB  * THE  NUMBER  OF  SUBNETWORKS  CREATED  IN  THE  CURRENT  STAGE 
TSNSUB  » THE  TEMPORARY  SUBNETWORK  BEING  USED  IN  STAGE 
SUBROUTINE 

TSUBN  » THE  NUMBER  OF  THE  SUBNETWORK  CURRENTLY  BEING  CONSIDERED 
TYPESN  ■ THE  TYPE  OF  SUBNETWORK  BEING  CONSIDERED! 

1 a BUNDLE  SUBNETWORK  2 > CUT  SUBNETWORK 

A MINIMUM  SUBNETWORK  IS  ONE  THAT  IS  NCT  COMPOSEO  OF  SMALLER 
SUBNETWORKS 

IMPLICIT  INTEGER*2  (A-ZI 

COMMON  NSUB«TNSUB,S«T,  SUBNET  , SOURC, SINKS, NARCSS 
DIMENSION  S<  300  I ,T( 300) , SUBNET!  300, 100) ,SCURC! 100) 

DIMENSION  SINKS!  100), NARCSS(iOO) 

THIS  PROGRAM  WILL  ACCCMODATE  A NETWORK  THAT  HAS  A MAXIMUM  OF 
300  NOOES,  A MAXIMUM  OF  300  ARCS  ANO  CAN  STORE  A MAXIMUM  OF 
100  SUBNETWORKS  IN  THE  BREAKUP  PROCESS 

NA  • THE  NUMBER  GF  ARCS  IN  THE  NETWORK  TO  BE  CONSIDERED 
NS  « THE  TOTAL  NUMBER  OF  POSSIBLE  SUBNETWORKS  IN  THE  BREAKUP 
PROCESS 

IN  THE  CASE  WHERE  NO  GOOO  ESTIMATE  CAN  BE  MAOE  OF  NS,  USE  NA 


*• 

3 

A 

5 

6 

7 

8 
9 

10 

LI 

12 

13 

14 

15 

16 
17 
IB 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
59 
59 


C FOR  ALL  DIMENSIONS 

C THE  ARRAY  DIMENSIONS  ARES 

C SINAI  .TINA), SUBNET! NA, NS) , SOURCI  NS) , SI NKS I NS ) , NARCSS INS) 

C 

C IN  SUBROUTINE  BUNDLE  THE  ARRAY  DIMENSIONS  ARE* 

C BDNUM(NS), CHECKINS) 

C 

C IN  SUBROUTINE  CUT  THE  ARRAY  DIMENSIONS  ARE: 

C ORIGIN!  NS) .POST  I NS) ,TARC I NS) ,RCUT INA) 

C 

C READ  THE  NETWORK  IN 

C 

6J00  CALL  NETIN  ILNODEN) 

STEP=1 
TYPE  SN3 1 
T SU3N* 1 

CALL  BUNDLE  ILNOOEN.TSUBN) 

C 

C PRINT  OUT  THIS  STAGE  0C  THE  BREAKUP 

C 

CALL  STAGE  I TYPESN.TSUBN) 

CTNSUB*TN$UB 
GO  TO  70 
C 

C FIND  THE  NUMBER  OF  NEWLY  CREATED  SUBNETWORKS 

C 

10  CTNSUR-0 

TSUBN^NSUa-LUOP 
TYPESN  = 1 
DO  20  1*1. LOOP 

C 

C FIND  THF  NEXT  SUBNETWORK  TO  BE  FURTHER  SUBDIVIDED 

C 

TSUBN*TSUBN«-l 

C 

C FIND  THE  LARGEST  NOOE  NUMBER  IN  SUBNETWORK  TSUBN 

C 

CALL  NQDER  I TLNDEN, TSUBN) 

C 

C FIND  THE  BUNDLE  SUBNETWORKS 

C 

CALL  BUNDLE  I TLNDEN, TSUBN) 

C 

C PRINT  OUT  THIS  STAGE  OF  THE  BREAKUP 

C 

CALL  STAGE  I TYPESN.TSUBN) 

C 

C IF  THERE  IS  ONLY  ONE  BUNDLE  FOUND  IN  SUBNETWORK  TSUBN,  WE  ARE 

C FINISHED.  PRINT  OUT  ITS  COMPONENT  ARCS 

C 

IF  ITNSUB.FQ.i)  CALL  ENOSNT  (TSUBN) 

C 

C COUNT  THE  NEW  NUMBER  OF  SUBNETWORKS  CREATED 

C 

20  C TNSUB*TNSUB+CTNSUB 
C 

C IF  ALL  SUBNETWORKS  ARE  IN  THEIR  SMALLEST  FORM,  WE  ARE  FINISHED 

C 

IF  ICTNSUB.EQ.O)  GO  TO  90 
GO  TO  70 


60 

61 

62 

63 

64 

65 

66 
b7 
63 

69 

70 

71 

72 

73 

74 

76 

77 

78 

79 

80  , 
81 
82 

83  * 

84  | 

85 

86 

87 

88 
39 

90 

91 

92 

93 

94 

95 


96 

97 

98 

99 
100 
101 
102 

103 

104 

105 


106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


80 


C 

C 

C 


c 

c 

c 


c 

c 

c 

c 

c 

c 

55 


c 

c 

c 


c 

c 

I C 

c 


c 

c 

c 

30 

c 

c 

c 


85 

70 

900 

C 

C 

c 

90 

100 

666 

9000 


L00P*CTNSU8 

CTNSUB*0 

TSUBN*NSUB-L00P 

PINO  THE  NUMBER  OF  NEWLY  CREATED  SUBNETWORKS 

TYPESN*2 
DO  30  1*1, LOOP 

FIND  THE  NEXT  SUBNETWORK  TO  BE  FURTHER  SUBDIVIDED 

T SU8N*TSU8N+ 1 
IF  (TSUBN.EQ.l)  GO  TO  55 

FIND  THE  LARGEST  NODE  NUMBER  IN  SUBNETWORK  TSUBN 

CALL  NODER  l TLNDEN, TSUBN) 

FIND  THE  CUT  SUBNETWORKS 

IF  (TSUBN.EQ.l)  TLNOEN*LNODE N 
CALL  CUT  I TSUBN, TLNDEN) 

PRINT  OUT  THIS  STAGE  OF  THE  BREAKUP 

CALL  STAGE  ( TYP E SN , TSUBN) 

'F  THFRE  ARE  NO  CUTS  FOUND  IN  SUBNETWORK  TSUBN,  WE  ARE  FINISHED 
PRINT  OUT  ITS  COMPONENT  ARCS 

IF  (TNSUB.EQ.l)  CALL  ENDSNT  (TSUBN) 

COUNT  THE  NEW  NUMBER  OF  SUBNETWORKS  CREATED 

CTNSUB*TNSUB+CTNSUB 

IF  ALL  SUBNETWORKS  ARE  IN  THEIR  SMALLEST  FORM,  WE  ARE  FINISHED 

IF  (CTNSUB.EQ.O)  GO  TO  90 
GO  TO  70 
LOOP*CTNSUB 
GO  TO  10 

WRITE  (6,900)  STEP 

FORMAT  ( 1H l , 5X, * STAGE' , 13,  • BREAKUP*) 

STEP*STEP*1 

LET’S  GO  BACK  TO  THE  APPROPRIATE  LOOP  FOR  THE  NEXT  STAGE 

GO  TO  (80,85) ,TYPESN 
CONTINUE 

READ  ( 5, 100,EN0*666) 

FORMAT  (13) 

GO  TO  6000 
CONTINUE 
WRITE  (6,9000) 

FORMAT  (1H1) 

RETURN 

ENO 


120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164  * 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 


SUBROUTINE  NETIN  ( LNODEN) 


178 


24 


•I 

t 


200 

210 

220 


C 

C 

C 

C 

C 

C 

C 

C 


\ 

i 


240 

C 

C 


C 

i 10 

; 


t 


ft 


IMPLICIT  INTEGERS  (A-Z) 

COMMON  N SUB,  TNSUB,S,T, SUBNET, SOURC,  S I NKS»NARCSS 
DIMENSION  S(  300) *T(  3001 1 SUBNET  I 300, 100) , SOURC I 100) 
DIMENSION  SINKS! 100) tNARC SSI  100) 


ZEROIZE  SOURC  ARRAY 


DO  20  I«l,50 
SOURC! I ) *0 

READ  IN  THE  INITIAL  NETWORK  LIMITS 

REAO  (5,100)  ARCS, SOURCE, SINK, LNODEN 
FORMAT  (413) 

WRITE  (6,200) 

FORMAT  (1H1,5X,* INPUT  STAGE*) 

WRITE  (6,210)  ARCS, SOURCE, SINK, LNQOEK 

FORMAT  ( 1H0, /6X, * THE  INITIAL  NETWORK  HAS* , 12X, 14, • ARCS', //6X, 

* • THE  SOURCE  IS  NODE  NUMBER* • 11X , I 3, / /6X ,* THE  SINK  IS  NODE  NUMBER', 
*13X, I3.//6X, • THE  LARGEST  NODE  IS  NODE  NUMBER* ,5X, 13) 

WRITE  (6,220) 

FORMAT  ( IHO, 5X, ' THE  INITIAL  NETWORK  AS  READ  IN  IS:',///6X, 

* ' ARC  NUMBER*, 5X, 'ORIGIN  NODE ' , 5X , • TERMI NAL  NODE') 

READ  IN  EACH  ARC  AND  ITS  STARTING  AND  TERMINATING  NODES 

THE  ARCS  AND  NODES  MAY  BE  NUMBERED  ANY  WAY  AND  READ  IN  IN  ANY 

ORDER 

I = ARC  NUMBER 

S « THE  NODE  NUMBER  FOR  THE  START  OF  AN  ARC 
T * THE  TERMINAL  NOOE  OF  AN  ARC 

DO  10  J*l, ARCS 

READ  (5,100)  I , S( I ) ,T ( I ) 

WRITE  (6,240)  I,S(I),T(!) 

FORMAT  (1H  ,8X,I3,13X,I3, 14X»I  3) 

CREATE  THE  FIRST  SUBNETWORK 

SUBNET  (J,1)*I 

SOURC! 1)«S0URCE 

SINKS! 1)»SINK 

NARCSS(n«ARCS 

NSUB«i 

RETURN 

END 

SUBROUTINE  NOOER  ( TLNDEN, TSUBN) 

FINDS  LARGEST  NODE  NUMBER  IN  THE  SUBNETWORK  TSUBN 
IMPLICIT  INTEGER*2  (A-Z) 

COMMON  NSIJB,  TNSUB  »S»T» SUBNET  , SOURC , SINKS,  NARCSS 
DIMENSION  S(  300) ,T( 300) , SUBNET! 300, 100) ,SOURCIiOO) 

DIMENSION  SINKS! 100), NARCSS(IOO) 

ARC S*NARCSS( TSUBN) 

TLNOEN-O 
DO  20  J»1,ARCS 
A»SU8NET(J, TSUBN) 

M*S( A) 

N*T (A) 


179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207  i 

208 


209 

210 
211 
212 

213 

214 

215 


216 

217 


218 

219 

220 
221 
222 
223 


224 

225 

226 
227 
220 

229 

230 

231 
2 32 

233 

234 

235 

236 

237 


MAXNO=N  238 

IF  (M.GT.N)  MAXND=M  239 

IF  (MAXNO.GT.TLNDEN)  TLNDEN*MAXND  240 

20  CONTINUE  241 

RETURN  242 

FNO  243 


SUBROUTINE  ENOSNT  ( TSU8N)  244 

C 245 

C PRINTS  SMALLEST  BREAKDCWN  OF  SUBNETWORK  TSUBN  246 

C 247 

IMPLICIT  INTEGER*2  (A-ZI  248 

COMMON  NSUB. TNSUB,S,T, SUBNET ,SOURC , S INKS  , NARCSS  249 

DIMENSION  S(300),T(300),SUBNET(300,100), SOURC 1 100)  250 

DIMENSION  SINKS!  100)  ,NARCSS(100)  251 

SOURCE*SOURC( TSUBN)  252 

SINK=SINKSITSUBN»  253 

WRITE  ( 6 f 100 ) TSUBN  254 

100  FORMAT  I 1H0,//16X,  ‘SUBNETWORK  *,13,'  IS  A MINIMUM  NETWORK 16X,  255 

* 1 IT  IS  COMPOSED  OF  : • 1 2 56 

WRITE  (6,200)  SOURCE, SINK  257 

200  FORMAT  ( 1H0 , 19X  ,•  SOURC  E NOOE  = ' , 13*  //20X,  'SINK  NODE  * M3)  258 

M*NARCSS(TSUBN)  259 

WRITE  (6,400)  260 

400  FORMAT  (IH0,19X,'ARC'«2X«'S( ARC ) ,,2X,,T(ARC)'I  261 

DO  10  1*1, M 262 

N*SU8NET(I,TSUBN)  263 

10  WRITE  (6,300)  N,S(N),T(N)  264 

30u  FORMAT  ( IH  ,19X,I3,3X,I3,5X,I3)  265 

C 266  | 

C SET  TNSUB-0  SO  THAT  THE  REMAINING  NUMBER  OF  SUBNETWORKS  DOESN'T  267 

C INCLUDE  THIS  MINIMUM  SUBNETWORK  263 

C 269 

TNSUB=0  270 

RETURN  271 

END  272 


C 

C 

C 


C 

C 

C 

C 

C 

c 


100 


300 


SUBROUTINE  STAGE  ( TYPESN, TSUBN) 

PRINTS  OUT  THE  CURRENT  STAGE  OF  BREAKUP 
IMPLICIT  INTEGERS  (A-ZI 

COMMON  NSUB,TNSUB,S,T, SUBNET, SOURC, SINKS, NARCSS 
DIMENSION  S( 3001  ,T(  300) , SUBNET (300, 100) , SOURC! 100) 

DIMENSION  SINKS!  100) , NARCSS! 100) 

TNSUB*THE  NUMBER  OF  NEW  SUBNETWORKS  RESULTING  FROM  THE  BREAKUP 
IN  THIS  STAGE 

TNSUB*l  IMPLIES  NO  BREAKUP  OCCURRED  IN  THIS  STAGE 

IF  (TNSUB.EQ.l)  GO  TO  600 
WRITE  (6,100)  TSUBN 

FORMAT  ( 1H0,/6X,  'SUBNETWORK  ',13,*  IS  COMPOSED  OF  SUBNETWORKS:  •) 

T SN  SUB*N SU8-TNSUB 

M«TSNSUB+i 

WRITE  (6,300)  ( I , IsM,NSUB) 

FORMAT  (1H0.5X, 20(13, «,')) 

IF  (TYPESN. EQ.il  GO  TO  60 
WRITE  (6,400) 


273 

274 

275 

276 

277 

278 

279  . 
230 
281 

2 82 

283 

284 

285 

286 

287 

288 
239 

290 

291 

292 

293 

294  * 

295 


o o o o oooo  oooo  ooo  ooo  o o o o 


90J  FORMAT  ( 1H0 , 5X,  • IN  SERIES')  29b 

GO  TO  600  297 

60  WRITE  (6,500)  298 

500  FORMAT  ( 1H0 , 5X , ' IN  PARALLEL')  299 

600  RETURN  300 

END  301 

SUBROUTINE  BUNDLE  ( LNODEN , TSU8N)  302 

303 

BUNDLE  IDENTIFIES  PARALLEL  SUBNETWORKS  CONNECTING  DESIGNATED  304 

SOURCE  AND  SINK  305 

306 

IMPLICIT  I NTEGER*2  « A-Z)  307 

COMMON  NSUB, TNSUB  »S  ?T • SUBNET ,S0URC , S I NKS  , NARCSS  308 

DIMENSION  S(300)«T(300), SUBNET  I 300, 100) , SOURC ( 100)  309 

DIMENSION  SINKS(IOO) ,NARCSS(100) , BDNUM(  300 ) , CHECK! 300)  310 

311 

GROUP  NODES  INTO  BUNOLES  312 

313 

SOURCE3 SOURC ( TSUBN)  314 

SINK*SINKS(TSUBN)  315 

NUMBD* 1 316 

317 

ZEROIZE  THE  BDNUM  ARRAY  318 

319 

DO  10  1= 1, LNODEN  320 

10  BDNUM ( I ) *0  321 

K*SUBNET(1 , TSUBN)  322 

M-S(K)  323 

N-T(K)  324 

BDNUM(M)«1  325 

BDNUM (N ) «1  326 

MAXND®N  327 

IF  (M.GT.N)  MAXND*M  328 

ARCS*NARCSS( TSUBN)  329 

IF  IARCS.EQ.1)  GO  TO  515  330 

DO  1 K=2» ARCS  331 

I«SUBNET(K,TSUBN)  332 

BDNUM  ( SOURCE  )-*0  333 

BDNUM! SINK >*0  334 

M=  S ( I ) 335 

N=T ( I ) 336 

IF(M.GT.MAXND)  MAXND-M  337 

IF(N.GT.MAXND)  MAXND=N  338 

339 

IF  AT  LEAST  l NODE  CN  THE  ARC  HAS  NOT  BEEN  ASSIGNED  TO  A BUNDLE  340 

GO  TO  2 341 

342 

I F( BDNUM ( M ) • EQ« 0 ) GO  TO  2 343 

344 

IF  ONLY  THE  TERMINAL  NOOE  ON  THE  ARC  HAS  NOT  BEEN  ASSIGNED  TO  A 345 

BUNDLE,  GO  TO  3 346 

347 

IF(BONUMIN).EQ.O)  GO  TO  3 348 

349 

IF  BOTH  NODES  ON  THE  ARC  HAVE  BEEN  ASSIGNED  TO  THE  SAME  BUNDLE  350 

EVERYTHING  IS  OKAY,  GO  TRY  ANOTHER  ARC  351 

352 

IF! BDNUM(N) •EQ«BDNUM( M)  ) GO  TO  1 353 

354 


C 


■-! 

t 


! 

c 

✓ 

27 

IF  THF  NODES  CN  THE  ARC  ARE  ASSIGNED  TO  DIFFERENT  BUNDLES, 

355 

1 

c 

THEN  THESE  TWO  BUNDLES  SHOULD  BE  POOLED 

356 

c 

357 

• 

I F( BDNUHIN) .LT. BDNUMI M J > GO  TO  6 

358 

c 

359 

c 

POOL  BUNDLES 

360 

c 

THE  BUNDLE  WITH  THE  LARGER  BUNDLE  NUMBER  IS  POOLED  INTO  THE 

361 

c 

BUNDLE  WITH  THE  SMALLER  BUNDLE  NUMBER 

362 

c 

THE  BUNDLE  NUMBERS  OF  ALL  BUNDLES  ARE  ALL  ADJUSTED 

363 

c 

364 

MAXBD*8DNUM|N I 

365 

MINBO«BONUM(M) 

366 

GO  TO  7 

367 

6 

M AXB0*BDNUM(  M ) 

368 

MINBD=BDNUM(N) 

36Q 

7 

00  5 J^ltMAXND 

370 

B«BDNUM< J) 

371 

IF  (B.EQ.MAXBD)  BDNUMt JJ=MI NBD 

372 

IF  (B.GT.MAX80)  BDNUMI J 1 "BDNyMI J 1 -1 

373 

5 

CONTINUE 

3 74 

NUMBD*NUMB0-1 

375 

GO  TO  1 

376 

. c 

377 

c 

IF  BOTH  NODES  ON  THE  ARC  ARE  UNASSIGNED,  GO  TO  4 WHERE  A NEW 

378 

c 

BUNOLE  IS  CREATED 

379 

c 

360 

2 

I F | BDNUMI N I . EQ.O)  GO  TO  4 

381 

c 

382 

c 

ASSIGN  THE  ORIGIN  NCDE  OF  THE  ARC  TO  THE  BUNDLE  CONTAINING  THE 

383 

c 

TERMINAL  NODE 

384 

c 

385 

BDNUMI M) "BDNUMI  N) 

386 

1 t: 

GO  TO  1 

387 

c 

388 

c 

ASSIGN  THE  TERMINAL  NODE  OF  THE  ARC  TO  THE  BUNDLE  CONTAINING 

389 

c 

THE  ORIGIN  NOOE  OF  THE  ARC 

390 

1: 

c 

391 

I- 

3 

BDNUMI N) "BDNUMI  Ml 

392 

j' 

- 

GO  TO  1 

393 

" j 

c 



394 

1 

c 

CREATE  A NEW  BUNDLE 

395 

| 

c 

396 

‘ 

4 

NUMBQ»NUMBD*1 

397 

1 

BDNUMI Ml »NUMBD 

398 

' 

BDNUMI NI*NUMBD 

399 

I 

1 

CONTINUE 

400 

515 

CONTINUE 

401 

i 1 

BONUMI SINK  I *0 

402 

c 

403 

c 

IF  WE  ONLY  HAVE  1 BUNDLE  FROM  THE  SUBNETWORK,  WE  ARE  FINISHED 

404 

I’1' 

c 

405 

IF  (NUMBD.EQ.l)  GO  TO  219 

406 

c 

407 

c 

ZEROIZE  CHECK  ARRAY 

408 

c 

409 

1 £' 

00  290  I "1 , ARCS 

410 

i ; 

290 

CHECK! I 1-0 

411 

I • 

L-0 

412 

c 

413 

c 

THE  NODES  ARE  IN  BUNDLES.  PUT  THE  ASSOCIATED  ARCS  INTO 

414 

C APPROPRIATE  PARALLEL  SUBNETWORKS  *15 

C *16 

DO  33  1 = 1, NUMBD  *17 

SUM ARC  *0  *18 

NSUB=NSUB* 1 *19 

00  3*  K=  I* ARCS  *20 

M=SUBNET(K , TSUBN ) *21 

N=S(M)  *22 

C *23 

C SOURCE  AND  SINK  HAVE  BUNDLE  NUMBER  0 *2* 

C *25 

IF  ( N .EO. SOURCE ) N=T( M)  *26 

IF  IBONUMINJ.EQ.II  GO  TO  239  *27 

IF  (N.EQ.SINKI  GO  TC  229  *28 

GO  TO  3*  *29 

C *30 

C SPECIAL  CASE : BUNOLE  HAS  ONLY  2 NODES:  SOURCEt  SINK.  *31 

C PUT  ALL  ARCS  THAT  ARE  PARALLEL  SUBNETWORKS  BY  THEMSELVFS  INTO  *32 

C THE  CHECK  ARRAY  *33 

C *3* 

229  DO  291  J *1 , K *35 

W=CHECK(  J)  *36 

IF  (W.EQ.O)  GO  TO  292  *37 

IF  (M.EQ.W)  GO  TO  3*  *38 

291  CONTINUE  *39 

.292  CHECK ( J ) =M  **0 

GO  TO  3*  **l 

C **2 

C THIS  ARC  IS  IN  THE  BUNDLE  It  HENCE  IT  IS  IN  THE  ITH  NEW  *43 

C SUBNETWORK  44* 

C 445 

239  SUMARC-SUMARC+1  446 

SUBNET (SUMARCtNSUB I 4*7 
3*  CONTINUE  *48 

C 449 

C CREATE  NEW  SUBNETWORKS  450 

C 451 

C 452 

C IF  THIS  BUNDLE  HAS  NO  NODES,  PUT  AN  ARC  SUBNETWORK  INTO  *53 

C SUBNET  (l.NSUB)  *5* 

C 455 

IF  (SUMARC.EQ.OI  GO  TO  333  *56 

343  NARCSS(NSU8I*SUMARC  *57 

SOURC INSUB) "SOURCE  458 

SINK SINS UB)*SINK  459 

GO  TO  33  460 

C *61 

C STORE  THE  SUBNETWORKS  THAT  HAVE  ONLY  SOURCE  AND  SINK  NODES  *62 

C 463 

333  SUMARC*1  464 

L=L*  l *65 

SUBNET(1,NSUB)>CHECK(L)  *66 

GO  TO  3*3  *67 

33  CONTINUE  *68 

*219  TNSUB-NUMBD  *69 

RETURN  *70 

ENO  *71 

SUBROUTINE  CUT  I TSUBN, LNODEN)  *72 

C *73 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


29 

CUT  IDENTIFIES  CUT  NOOES  EXCLUDING  THE  DESIGNATED  SOURCE  ANO 
SINK 

CUT  ALSO  IDENTIFIES  THE  CUT  GROUPS!  THAT  IS,  THE  SUBNETWORKS 
WHICH  ARE  IN  SERIES  ANO  CCNNECTED  BY  THE  CUT  NODES 

IMPLICIT  INTEGER*2  (A-Z» 

COMMON  NSUB, TNSUB,S,T, SUBNET ,SUURC, SINKS, NARCSS 

DIMENSION  S ( 300 ) »T(  300 ), SUBNET  I 300, 100) , SOURC  ( 100 i , TARC ( 300 ) 

DIMENSION  SINKS! 100) ,NARCSS< 100 ) ,ORI GIN! 300) , POST ( 300) , RCUT ( 100) 

FIND  THE  CUT  NODES 

NCUT  IS  THE  NUMBER  OF  CUT  NODES  FOUND  THUS  FAR 
NCUT  *0 

ARC  S=NARCSS l TSUBN) 

SOURCE=SDURC 1 TSUBN) 

SINK*SINKS( TSUBN) 

THE  DO  LOUP  DOWN  TO  STATEMENT  NUMBER  1 DETERMINES  THE  CUT  NODES 
DO  1 K»l, LNODEN 

CHECK  TO  SEE  IF  NODE  K IS  ACTUALLY  IN  THE  SUBNETWORK 

DO  20  J*  1 , ARC  S 
2-SUBNETI J, TSUBN) 

IF(Sm.EQ.K)  GO  TO  21 
IF  (T(Z).EQ.K)  GO  TO  21 
CONTINUE 

NODE  K IS  NOT  IN  THIS  SUBNETWORK 

GO  TO  1 
CONTINUE 

NODE  K IS  IN  THIS  SUBNETWORK 

IFIK.EQ. SOURCE)  GO  TO  1 

IFIK.EQ. SINK)  GO  TO  l 

NTARC-0 

DO  2 J*1 , ARCS 

Z*  SUBNET!.),  TSUBN! 

IF  ( S(Z).EQ.K)  GO  TO  2 
IF  IT(Z).EQ.K)  GO  TO  2 
NTARC-NTARCH 

TARC (N TARC) “SUBNET!  J, T SUBN) 

CONTINUE 

TARC  IS  THE  SUBNETWORK  WITHOUT  THE  ARCS  INVOLVING  NODE  K 
IF  TARC  CONTAINS  A PATH  FROM  THE  SOURCE  TO  THE  SINK*  THEN  NODE 
K IS  NOT  A CUT  NODE 
OTHERWISE,  K IS  A CUT  NODE 

OR  I G IN!  I )» SOURCE 
NORIG* 1 
CONTINUE 
NPOST-O 

IF  THERE  ARE  NO  ARCS  IN  THE  TARC  ARRAY,  K IS  A CUT  NODE 


474 

475 

476 

477 
4 78 
4 79 

480 

481 

482 
433 

484 

485 

486 

487 

488 

3 r 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501  | 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 


30 


IF  (NTARC.EQ.OJ  GO  TO  44 

FIND  ALL  NOOES  WHICH  CCME  AFTER  AN  ORIGIN;  PUT  THEM  IN  POST 
DO  4 1*1, NORIG 
DO  5 J=l,NTARC 
Y=OR I G IN ( I ) 

Z*TARC ( J ) 

U=S(Z) 

V-T I Z ) 

IF  ( U.NE.Y ) GO  TO  5 
IF  (V.EQ.SINK)  GO  TO  1 


IF  WE'VE  REACHED  THE  SINK,  NODE  K IS  NOT  A CUT  NOQE 


IF(NPOST.GE.l)  GO  TO  8 

NPOST=NPOSm 

POST (NPOST )*V 

GO  TO  5 

DO  9 L*i, NPOST 

IF  THIS  TERMINAL  NODE  IS  ALREADY  A POST,  LET'S  IGNORE  IT 

IF  (POSTIL I.EQ.V)  GO  TO  5 

CONTINUE 

NPOST-NPOSm 

POST (NPOST )*V 

CONTINUE 

CONTINUE 

IF  THERE  ARE  NOW  NO  POSTS,  NODE  K IS  A CUT  NODE 

IF(NPOST.NE.O)  GO  TO  13 

NCUT-NCUTM 

RCUT(NCUT)»K 

GO  TO  l 

NORIG-NPOST 

DO  14  L«l,NORIG 

THESE  ARE  NOW  OUR  NEW  ORIGINS 

ORIG  IN ( L I* POSTIL  I 

CHECK  THE  NEW  ORIGINS  FOR  THEIR  POSTS 

GO  TO  11 
CONTINUE 

IF  (NCUT.EQ.O)  GO  TO  32 
NSUB*N SUB+1 
SOURC(NSUB)-SOURCE 
OR  I GIN  ID- SOURCE 

NOW  WE  NEED  TO  FIND  THE  COMPONENTS  OF  THE  SERIES  SUBNETWORKS 
THAT  ARE  SEPARATED  BY  THE  CUT  NODES 

SUMARCO 
NORIG* i 
NPOST-O 


534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 


31 


I 


C 

C ZEROIZE  POST  ARRAY 

C 

on  55  1=1*  LNOOEN 
55  POST ( I )*0 

DO  24  1*1, NORIG 
YsORIGINlI) 

00  25  J=l, ARCS 
Z=SUBNETI J « TSUBN) 

C 

C ALL  ARCS  BEGINNING  AT  THIS  ORIGIN  GO  INTO  THE  NEW  SUBNETWORK 

C 

IE  (S(Z).NE.YI  GO  TO  25 
SUM ARC* SUM ARC *1 
SUBNET  ( SUMARC,NSUB!=Z 
C 

C TIZ)  WILL  BE  A NEW  ORIGIN  IF  IT  ISN'T  A REPEAT  OF  A CURRENT 

C ORIGIN 

C 

C CHECK  TO  SEE  IF  IT  IS  A REPEAT 

C 

00  30  K*l, NORIG 
X*  T I Z ) 

C 

C IF  TIZ)  IS  A REPEAT  OF  A CURRENT  ORIGIN,  LET'S  IGNORE  IT 

C 

IFIX.EO.ORIGINIK))  GO  TO  25 
30  CONTINUE 

NPOST=NPOST+i 
POST  INPOST ) *X 

C 

C CHECK  TO  SEE  IF  ANY  POST  IS  A CUT  NODE 

C IF  IT  IS,  REPLACE  IT  WITH  A ZERO 

C 

DO  52  W* l,NCUT 
0*RCUT ( W) 

IF  IX.NE.O)  GO  TO  52 
POSTINPOST  1*0 
NP0ST*NP0ST-1 

C 

C THIS  CUT  NODE  IS  THE  SINK  OF  THE  SUBNETWORK  UNOER 

■C  CONSIDERATION  ANO  THE  SOURCE  OF  THE  NEXT  SUBNETWORK  TO  BE 

C CONSIDERED. 

C 

SINKSINSUBI-O 

SCURCINSUB*ll«D 

C 

C CHECK  TO  SEE  IF  ANY  POSTS  ARE  REPEATED  IN  THE  POST  ARRAY.  IF 

C THEY  ARE,  REOUCE  THE  NUMBER  OF  POSTS  TO  WHERE  THERE  ARE  NO 

C REPEATS 

C 

52  CONTINUE 
C 

C IF  WE  HAVE  ONE  OR  LESS  POSTS,  THERE  ARE  NO  ADJUSTMENTS  OF  THE 

*C  POST  ARRAY  TO  BE  MADE;  LET'S  CONTINUE 

C 

53  IF  INPOST.LE.il  GO  TO  25 
POSTCK-NPOST-l 

DO  7 K»I,POSTCK 

IF  I POST  IK ) .E0« POST ( NPOST) ) GO  TO  61 
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651 
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32 


61 

25 

24 


8 


31 

32 


CONTINUE 

G(J  TO  25 

POSTINPOSTJ-O 

NPCST=NP0ST-1 

CONTINUE 

CONTINUE 

IF  WE  HAVE  NO  POSTS  LEFT#  WE 


HAVE  FOUND  AIL  OF  THIS  SUBNETWORK 


IF  (NPOST.EQ.OI  GO  TO  34 
NORIG* NPOST 
00  28  L*l, NORIG 
ORIGINIL)* POSTIL  I 
GO  TO  23 

NAKCSSINSU8) * SUM ARC 
NSUB=NSUBF 1 
X*  SOURC 1 NSUB ) 


IF  THE  SOURCE  OF  NSUB  IS  NOT  A CUT  NODE,  WE  NEED  TO  ADJUST  NSUB 
AND  GO  BACK  TO  THE  MAIN  PROGRAM  FOR  THE  NEXT  STAGE  OF  THE 
BREAKUP 


IF  (X.EQ.O)  GO  TO  31 

ORIGIN! 1 1*  X 

GO  TO  39 

N $UB*N  SUB- 1 

SINK  SINS UB)*SINK 

TNSU8*NCUm 

RETURN 

END 


654  3 

655  1 

656 

657 
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669 

670  ' 

671 


672 
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674 
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677 

678 

679 

680 
681 
682 
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This  short  program  will  determine  whether  a given  network  is  acyclic 
or  contains  loops  (cycles) . The  program  examines  each  node  and  Indicates 
whether  or  not  the  node  is  part  of  a loop.  If  a node  is  part  of  a loop, 
the  number  of  activities  in  the  loop  is  also  indicated. 

The  basic  steps  in  determining  whether  or  not  the  INODE-th  node  is 
part  of  a loop  are  as  follows: 

(1)  Identify  all  activitiea  whose  termanal  node  is  the  INODE-th 
node.  Let  A be  the  set  of  all  origin  nodes  for  these  activities. 

(2)  If  INODE  is  in  A,  tha  INODE-th  node  is  part  of  a loop  and  stop. 

(3)  Identify  all  activitiea  whose  terminal  node  is  In  A.  Redefine 
A to  be  the  sat  of  origin  nodes  for  these  activities.  If  A is 
now  empty,  stop  and  the  INODE-th  node  is  not  part  of  a loop. 

If  A is  not  empty,  return  to  step  2.  ' 

Specific  Input  Instructions: 

Card  1.  Col.  1-3:  The  number  of  activities  in  the  network.  Format  13). 

Col.  4-6:  The  largest  node  number  in  the  network.  Format  (13). 

For  each  activity  one  card  with: 

Col.  1-3:  The  activity's  origin  node  number,  Format  (13). 

Col.  4-6:  The  activity's  terminal  node  number,  Format  (13). 

•i 

The  activities  may  be  inputted  in  any  order.  The  nodes  may  be  num- 
bered in  any  manner;  however,  the  program  is  sure  efficient  if  the  set  of 
node  numbers  contains  only  the  numbers  1 through  M where  N is  the  number 
of  nodes  in  tha  network. 


Dimension  Restrictions 


This  program  is  written  in  FORTRAN  G.  The  current  dimensions  will 
allow  a network  with  300  activities  and  200  nodes  to  be  considered. 


SAMPLE  NETWORK 


Node  1 is  not  looped 

Node  2 ie  looped.  There  ere  3 activities  in  the  loop. 
Node  3 is  looped.  There  are  3 activities  in  the  loop. 
Node  4 is  looped.  There  are  3 activities  in  the  loop. 


oooooo  ooooo 


PROGRAM  LISTING 
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PROGRAM  LOOP  1 

2 

THIS  PROGRAM  DETERMINES  WHETHER  THE  INPUTTED  NETWORK  IS  ACYCLIC  3 

OR  CONTAINS  NODES  WHICH  ARE  PART  OF  LOOPS  ICYCLESI  A 

5 

IMPLICIT  I NTEGER62  IA-Z)  6 

DIMENSION  A ( 200), 81  200) , HE ADI  300 ), TAIL ( 300 ) 7 

8 

THE  ARRAY  DIMENSIONS  ARE:  AIN ) , BIN) ,HEAD( M ), TA IL IM ) 9 

WHERE  10 

N = THE  NUMBER  OF  NODES  IN  THE  NETWORK  U 

M = THE  NUMBER  OF  ACTIVITIES  IN  THE  NETWORK  12 

13 

( L0QP*0  16 

READ  (5,100)  M, NMM  15 

100  FORMAT  (213)  16 

RE  ADI  5,101)  (TAIL!  I), HE AD ( 1),1*1,M)  17 

Ul  FORMAT  (213)  18 

WRI TE I 6, 2001 ) 19 

2001  FORMAT ( IH 1 ) 20 

WRI TE ( 6,2000)  ( I , TAI L ( 1 1 , HE ADI  I ) , I * 1, M)  21 


2000  FORMAT (•  ACTIVITY  TAIL  HEAD' / , (6X,  1 3, 7X, 1 3,6X , 1 3) ) 22 


WRITE  (.6,2001)  23 

C FORM  THE  1ST  HIERARCHY  26 

INODE  * 0 25 

80  HI ER  - 2 26 

INODE  * INODE  ♦ l 27 

J * 0 28 

00  l I-  1,M  29 

IF(HEAO(I).NE. INODE)  GO  TO  1 30 

J « J*1  31 

A(J)=TAIL(I)  32 

IF  (TAIL!  I). EQ. INODE)  GO  TO  998  33 

1 CONTINUE  36 

IF  IJ.EQ.O)  GO  TO  997  35 

IA*J  36 

J*0  37 

C FORM  THE  SUBSEQUENT  HIERARCHIES  38 

102  CONTINUE  39 

DO  2 11=1, IA  60 

DO  3 1*1, M 61 

IF  (HEAD!  I).NE.A(ID)  GO  TO  3 62 

IF(  TAIL!  I). EQ. INODE)  GO  TO  998  63 

IF(J.EQ.O)  GO  TO  60  66 

DO  10  K*l,J  65 

IF(  TAIL(t).EQ.BIK) ) GO  TO  11  66 

10  CONTINUE  67 

60  CONTINUE  68 

J*J+l  69 

81 J)*  TAIL!  I ) 50 

11  CONTINUE  51 

3 CONTINUE  52 

2 CONTINUE  53 

IF  (J.EQ.O)  GO  TO  997  56 

HIER  « HI ER* 1 55 

1 A* J 56 

J = 0 57 

00  20  1*1, IA  58 

20  A( I ) =B(  I ) 59 

“ "■‘''"•"'••““i.  - >■  aA.,- .wvv; i-i. 
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997 

2002 


1000 


GO  Tl  102 
CONTINUE 

WRITE! 6*2002 1 INODE 

FORMAT! • NODE  * « 15 1 * IS  NOT  LOOPED*  ) 

IF  ! INODE. NE.NMM  } GO  TO  80 
IF!  ILOOP.EQ. U GO  TO  50 
WR l TE I 6» 1000 ) 

FORMAT  (•  THERE  ARE  NO  LOOPS  IN  THIS  NETWORK  * I 


GO  TO 

50  MRI TE ( 6*  51 ) 

51  FORMAT! • THERE  ARE  NO  OTHER  LOOPS  IN  THIS  NETWORK*) 

GO  TO  999 

998  WRITE  (6.1001)  INODE, HIER 

1001  FORMAT  (•  NODE  • #15,*  IS  LOOPED.  THERE  ARE  *,I3,*  ACTIVITIES  IN  T 
♦HE  LOOP.* ) 

ILOOP*l 
INDDD«INODEM 
IF!  INDDO. NE.NMM)  GC  TO  80 
999  CONTINUE 

WR I TE ( 6.2001 ) 

STOP 

END 
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